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Abstract. 

, We investigate how the geometrical and physical properties of soft X-ray flare kernels 

change with their altitude above the photosphere. We analyze limb flares well observed by 
Yohkoh/S^T showing clear geometry with well separated loop-top kernels. Our analysis 
concerns relations between kernel size, plasma pressure, energy release and the kernel 
altitude. We define scaling laws describing how the sizes and its physical properties of 
kernels vary with the altitude above the photosphere. We interpret the observed relations 
in terms of the general magnetic structure of active regions. 
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, 1. Introduction 

The loop structures making up soft X-ray (SXR) flares generally show bright 
!>• \ regions at their tops, which we here call kernels. They form before maximum 

and in long-duration flares may last for hours, with new kernels forming 
\^ \ after the flare maximum. A kernel's surface brightness may be ten times as 

' bright as the rest of the flare loop. High-resolution images show the thermal 

. structure to be uniform (Feldman et al., 1995; Jakimiec et al., 1998) and 

\ about half the total flare emission measure is contained with the kernel. 

■ The classic model describes solar flare as a hydrodynamic reaction of 

dense chromospheric plasma to a sudden release of energy in the corona. 
This energy conducted down along the magnetic lines and transfered by 
^ . energetic electrons heats up the cool, dense chromosphere and forces it to fill 

magnetic flux-tube, where the energy release occurred (Bentley et al., 1994; 
Tomczak, 1997). In a case of simple magnetic structure the plasma should 
fill it uniformly, without local condensation. The observed kernels can be 
brighter due to higher density or temperature, and in a simple magnetic 
loop any kernel should disappear in a very short time due to expansion or 
conduction. Some kind of restriction efficiently preventing outflow of mass 
and energy even for hours must be present at the boundary of a kernel 
(Vorpahl, Tandberg-Hanssen and Smith, 1977). 

Loop-top kernels are recognizable in the the flare images from Skylab. 
Although widely commented on since the Yohkoh Soft X-ray Telescope 
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(SXT) observations (e.g. Acton et al., 1992; Doschek, Strong and Tsuneta, 
1995; Jakimiec et al., 1998), they are still not well understood. Recently 
proposed mechanisms are based on MHD turbulence (Jakimiec, 2002), fast- 
mode MHD shock (Yokoyama and Shibata, 1998) or magnetic trap at the 
top of the cusp structure (Karlicky and Barta, 2006). 

The ratio of the kernel soft X-ray emission to the total emission of the 
flare is in the range 0.35 to 0.45 for different flares, and moreover remains 
nearly constant for the whole period of kernel existence (Doschek, Strong and 
Tsuneta, 1995). This property makes the kernels a usefull tool to analyze the 
evolution of the whole flare. In this paper we focus on how the geometrical 
and physical properties of loop-top kernels depend on their altitude above 
the photosphere. To this work we utilize the observations of Yohkoh SXT 
instrument (Tsuneta et al., 1991), which form a large database of soft X-ray 
flares over the period from October 1991 to December 2001 (maximum of 
Cycle 22 to maximum Cycle 23). Data from the SXT allow us to determine at 
the same time both geometrical properties of flare kernels as well as physical 
properties of the plasma. 



2. Sample description 

To estimate a kernel altitude (h) we need to know its position on the SXT 
image and assume above which point of the photosphere it is situated (sub- 
kernel point). The Ha database and 'Masuda method' (Masuda, 1994) are 
the most helpful in locating this point. The uncertainty of its position is 
the main source of error in estimating the kernel altitude. The uncertainty 
is least at the limb, so we chosed to analyze events close to the limb. We 
searched the GOES database for flares with location within 10 degrees from 
the limb. Looking at SXT images of these flares we chose events with simple 
structure showing a single loop-top kernel or in the case of multi-kernel 
flares we chose the kernels well separated from other sources. Wc enlarged 
the sample with some flares at longitudes slightly grater than 10 degrees from 
the limb or located behind the limb. To obtain reliable physical properties 
of the emitting plasma we selected the flares observed with the use of the 
SXT thick aluminium (A112) and beryllium filters (Bel 19) in full angular 
resolution (2."45) during the flare maximum. 

As each flare evolves greatly with time, we compared the properties of 
flare kernels at the time of maximum emission in the SXT All 2 filter. We 
found 48 suitable events (see Table I) during the entire period of Yohkoh 
operations (Oct 1991 to Dec 2001). The sample consists of 8 X-class flares, 
24 M-class and 16 C-class events. 34 of them were located within 10 degrees 
from the limb, while for 6 flares the SXR footpoints seen on SXT images 
were located between 60 and 80 degrees of heliographical longitude. To this 
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Table I. The sample of selected flares. Remarks in the last column: A - flare with one kernel, B - flare 
with more than one kernel but only one analyzed, C - flare with two kernels, both analyzed. 



No 


date 


maximum 


GOES 
UT] X- 


ray class 


co-ordinates 


location 

NOAA active region 


remarks 


1 


17-NOV-91 


07:16 




Ml.l 


E87 


N13 




6929 


A 


2 


19-Nov-91 


09:32 




C8.5 


W64 


S13 




6919 


A 


3 


02-Dec-91 


05:01 




M3.6 


E87 


N16 




6952 


B 


4 


09-Dec-91 


09:44 




M4.1 


E89 


S05 




6966 


A 


5 


nQ-ner-QI 


23:49 




Ml.O 


E81 


S07 




6966 


A 


6 


13-Jan-92 


17:34 




M2.0 


W90 


S16 




6994 


A 


7 


30- Tan-92 


17:15 




Ml. 6 


E83 


S13 




7042 


A 


8 


06-Feb-92 


03:29 




M7.6 


W85 


N06 




7030 


c 


9 


17-Feb-92 


15:46 




Ml. 9 


W80 


N15 




7050 


A 


10 


28-.Tun-92 


05:14 




XI. 8 


W102 


N13 




7205 


B 


11 


28- liin-Q2 


14:24 




Ml. 6 


E90 


N14 




7216 


B 


12 


05-Jul-92 


12:04 




C4.1 


E82 


Sll 




7220 


A 


13 


12-Oct-92 


21:53 




C2.5 


W83 


819 




7303 


A 


14 


n2-Nov-Q2 


03:08 




X9.0 


W99 


S24 




7321 


A 


15 


21-Nov-92 


07:13 




C5.0 


W81 


816 




7341 


A 


16 


23-Nov-92 


13:59 




C4.0 


W84 


808 




7342 


A 


17 


29-Nov-92 


08:58 




C9.1 


W90 


827 




7345 


A 


18 


1 7-Feb-Q3 

J. 1 J. ^VJ iJKi 


10:40 




M5.8 


W88 


807 




7420 


A 


19 


02-Mar-93 


15:10 




C5.0 


E82 


S04 




7440 


A 


20 


15-Mar-93 


21:35 




M2.9 


W94 


802 




7440 


A 


21 


12-,Jun-93 


09:06 




C3.5 


W95 


811 




7518 


A 


22 


25-Jun-93 


03:22 




M5.1 


E84 


S09 




7530 


A 


23 


27-ScD-93 


12:12 




Ml. 8 


E86 


NIO 




7590 


A 


24 


28-Jan-94 


17:04 




Ml. 8 


W87 


N06 




7654 


A 


25 


27-Feb-Q4 


09:20 




M2.8 


W97 


NIO 




7671 


c 


26 


30- A US'- 94 


19:54 




C6.2 


E77 


809 




7773 


A 


27 


21-Apr-95 


13:41 




C5.1 


W75 


801 




7863 


A 


28 


17-SeD-97 


11:43 




Ml. 7 


W81 


N21 




8084 


A 


29 


2fi-Nov-97 


04:47 




C4.7 


E86 


N20 




8113 


A 


30 


08-May-98 


02:04 




M3.1 


W88 


816 




8210 


A 


31 


18-Aue'-98 


08:24 




X2.8 


E91 


N33 




8307 


A 


32 


18-Aue-98 


22:19 




X4.9 


E86 


N33 




8307 


A 


33 


1 Q_ Aii£i'-98 


14:26 




M3.0 


E79 


N33 




8307 


A 


34 


22-Nov-98 


06:42 




X3.7 


W76 


828 




8384 


A 


35 


22-Nov-98 


16:23 




X2.5 


W81 


828 




8384 


A 


36 


24-Nov-98 


02:20 




Xl.O 


W98 


829 




8384 


A 


37 


23-Dor-9S 


06:59 




M2.3 


E90 


N23 




8421 


B 


38 


25-Jul-99 


13:38 




M2.4 


W82 


N39 




8639 


B 


39 


21-Sep-99 


10:47 




C6.4 


W96 


825 




8692 


A 


40 


23-May-OO 


17:54 




C4.3 


W75 


821 




8996 


B 


41 


01-Jun-OO 


06:17 




M2.5 


E80 


N21 




9026 


A 


42 


30-Sep-OO 


23:21 




X1.2 


W93 


N08 




9169 


A 


43 


14-Oct-OO 


08:40 




Ml.l 


W81 


N02 




9182 


A 


44 


14-Oct-OO 


12:05 




C8.4 


W80 


N02 




9182 


A 


45 


28-Oct-OO 


07:10 




C9.7 


E80 


N08 




9212 


A 


46 


14-Nov-OO 


16:34 




Ml.O 


E82 


N13 




9233 


A 


47 


11-Mar-Ol 


08:56 








x^^^/02/2008; 


mt; p.3 


A 


48 


29-Oct-Ol 


01:59 




M1.3 


W88 


N13 






A 
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number we added 8 flares located no more than 12 degrees behind the 
hmb for which the loop-top kernels were not occulted by the limb. In our 
sample only two flares had two well-scparatcd loop-top kernels, allowing us 
to analyze both kernels. 6 flares had more than one kernel but only one of 
them could be analyzed. The rest of our events (i.e. 40) had only one kernel. 
Thus, a total of 50 flare kernels were analyzed. 



3. Analysis 

The geometrical parameters of interest here arc the size and altitude of 
each kernel, and the physical ones are plasma density, pressure and energy 
release rate in the kernel. The kernel size we defined by a 50% isophote (/50) 
relative to the brightest pixel in the A112 image. We measured the area. A, 
of the kernel projected on the plane of the sky. The kernel altitude, h, was 
calculated from the position of its centroid and coordinates of the sub-kernel 
point. The kernel centroid was estimated as a center of gravity in the sense 
of intensity distribution within the isophote /50. From the kernel area we 
summed the A112 and Bell9 flux to estimate the emission measure, EM, 
and temperature, T, of emitting plasma using the filter ratio method as 
described by Hara et al. (1992). 

The estimation of kernel volume is not straightforward. Usually the kernel 
image has a circular or elliptical shape. If we assume that its volume is 
ellipsoidal we can calculate it as F = (4/3)7ra6c where a, b and c are the 
three semi- axes. In general we can estimate two semi- axes from the kernel 
image, but the third dimension along the line-of-sight must be reasonably 
guessed. Generally the smaller of the two taken from the image is used to 
estimate the volume. This assumes that the main axis of ellipsoid is located 
in the plane of image. However, wc do not know definitely the orientation 
of the kernel relative to the plane of the sky, so the effective depth of the 
kernel can be greater. To avoid this underestimation we assume that the 
depth is close to the geometrical mean of the two axes seen on the image. 
This assumption gives the estimation of kernel volume as 

3 0r ' 

where A is the plane-of-sky kernel area. 

Having estimates of the kernel volume we can calculate the plasma mean 
density and pressure within the kernel as A^e = \jEMjV and Pe = 2A;sriVe 
where fcs is Boltzmann constant. One should remember that such estima- 
tions are sensitive to the influence of filling factor and in general should be 
treated as lower estimates. 
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The kernel area is defined by half the intensity of its brightest pixel, so 
the uncertainty in this pixel's flux directly affects the uncertainty A/50 of the 
isophote /50 and the uncertainty AA of the kernel area. Estimation of is 
also affected by the data pixelation. For small kernels the pixelation makes 
it impossible to estimate AA, because there may be no change of A within 
the range (I50 — AJ50, /50 + AJ50). To overcome this limitation we closely 
analyzed for each kernel how the number of pixels brighter than a given 
level changes with an isophote I in the range (/50 — 5 A/50, -^50 + 5 A/50). 
The step function obtained we fitted by a low degree polynomial. Using 
this polynomial we interpolated the area A at the /50 level and estimated 
the error A A for A/50. In the same way we calculated the flux emitted 
from within the isophote /50 in both SXT filters and an error of the flux. 
Our method allowed us to estimate the size of a kernel and its flux with 
precision greater than one pixel. 

The error AA affects not only the estimated volume but also all physical 
properties of kernel plasma: the total brightness in both fllters, and therefore 
the kernel temperature and emission measure, electron density and pressure. 
The error of kernel altitude comes from two uncertainties, one related to the 
kernel centroid position and second to the location of sub-kernel point. The 
second uncertainty is usually several times larger than the flrst. All errors 
mentioned affect the calculation of heating rate within the kernel. 



4. Results 

The flare kernel altitudes in our sample vary from 3.2 x 10* to 5.8 x 10^ cm. 
The distribution of altitudes decreases with rising altitude. This distribution 
is shown in Fig. 1 and may be described as exponential, N oc exp(—h/ho), 
where ho = (2.2 ± 0.6) x 10^ cm. 

Figure 2 presents the relation of the kernel size and altitude during the 
course of a few arbitrarily selected flares. As can be seen, individual flares 
may evolve in this diagram in very different ways, showing both correlation 
and anti-correlation as well as lack of any relationship. However, when we 
compare all the flares in our sample the general tendency becomes clear: 
the higher the kernel, the larger its size. Figure 3 illustrates this effect in 
an example of three arbitrary selected flares. The changes of kernel size 
and altitude during the course of all analyzed flares are shown in Fig. 4. 
To correctly compare all kernels we had to choose one specific moment in 
their evolution, which we chose the moment of kernel brightness maximum. 
Figure 5 presents the result of our analysis. The projected area of a loop-top 
kernel distinctly rises with its altitude. The relation is a power-law, Aoc h"', 
where the index n = 1.13 ± 0.04. This relation has power index 0.56 for the 
kernel mean radius, and 1.69 for the volume. 
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1x10^ 2x10^ 3x10^ 4x10^ 5x10^ 6x10^ 7x10^ 

h [cm] 

Figure 1. The observed distribution of loop-top kernel altitudes. 
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Figure 2. The changes of the projected kernel area and its altitude above the photosphere 
for a set of arbitrary selected flares. Individual kernels show different correlations but the 
overall tendency is visible. The numbers denote the flares listed in Table I. 
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SXT/AI12 17-Feb-1992 (9) SXT/AI12 06-Feb-1992 (8) SXT/AI12 23-Dec-1998 (57) 




860 880 900 920 940 960 980 920 940 960 980 1000 1020 1040 -1000-980 -960 -940 -920 -900 -880 

X (orcsecs) X (orcsecs) X (orcsecs) 



Figure 3. Three examples of flares with loop-top kernels at different altitudes. The ex- 
tension of each kernel is marked by a 50% isophoto relative to the brightest pixel in SXT 
A1.12 image. Each image covers the area of 140" x 140". Note the increase of the kernel 
size when going to higher events. 




10^ 10^ 10^° 

h [cm] 

Figure 4- The changes of the projected kernel area and its altitude above the photosphere. 
Measured kernel areas and altitudes at all times during all the flares analyzed are plotted 
here. The overall tendency is clearly visible, higher kernels are usually bigger. 
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Figure 5. Relation between the projected kernel area and its altitude above the photo- 
sphere for the time of maximum kernel brightness in the SXT A112 filter. The line fitted 
is the power-law with slope 1.13 ± 0.04. 

Our sample shows no relation between temperature or emission measure 
and the altitude of the kernel. Nor is there correlation with kernel area. 
Kernel emission measure correlates well with GOES X-ray class. Its median 
value rises from 3.1 X 10^^ cm~^ for C-class flares through 1.2 x 10^^ cm^'^ for 
M-class flares to 1.2 x lO^'' cm~^ in X-class events. The dependence of kernel 
temperature with GOES class is weaker but also evident in our sample. This 
relation has been often reported and discussed in the literature in similar 
form as the T—EM relation. In our sample it can be well described by power- 
law relation T oc EM^'^'^^^'^'^ (see Fig. 6). This dependence is weaker than 
others mentioned in literature. The reported power-law slopes vary from 
0.13 (Shibata and Yokoyama, 2002) to 0.23 (equivalent of the exponential 
relation in Feldman et al. (1996)). This weak dependence in our sample 
is because of the small sensitivity of SXT to hot plasma (Jakimiec et al., 
1998). Reconstructions of differential emission measure for flares generally 
show the presence of two components. One component contains the relatively 
cool plasma with temperatures between 5 MK and 10 MK and the other 
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10^ 10^ 10^ 10^ 

EM [10^^ cm"^] 

Figure 6. Observed relation between the mean kernel temperature and emission measure. 
The fitted line has slope 0.09 ± 0.02. Here and in following figures we apply different 
symbols for flares of different GOES class. 



component contains the hotter plasma with T between 15 MK and 25 MK 
(Kgpa ct al., 2005). In the presence of the cooler component, the hotter 
plasma only weakly contributes to the SXT signal. 

Figure 7 shows the relation between kernel density and altitude. Taking 
the whole sample we can hardly see any relationship, but if we divide the 
sample into the GOES classes the relation becomes more obvious. In a given 
subclass the kernel density drops with rising altitude. The observed power- 
law index of this relationship for C, M and X-class flares is equal — 1.12ib0.11, 
—0.82 lb 0.04 and —0.65 it 0.12 respectively. The difference between slopes for 
M and X-class kernels is barely significant. The kernels of C-class flares seem 
to follow a significantly steeper slope, but the reliability of its estimation 
is weak. This sub-sample is ill-spaced with height, as it includes only 3 
kernels with altitude higher than 2 x 10^ cm. For such a sample 'bootstrap' 
methods give more reliable estimation of slope. Using a bootstrap resampling 
method (Press, Teukolsky and Vatterling, 1992) we receive the same slope 
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but with error 0.60, which makes differences between all 3 slopes statistically 
unimportant. Similar relations are observed between plasma pressure and 
kernel altitude. In Fig. 8 we observe very similar power-law slopes for each 
GOES-class sub-sample (C: -1.07±0.12, M: -0.85±0.05, X: -0.69±0.14). 

These slopes seem to be defined by the observed geometrical properties 
of flare kernels. We do not observe any correlation between temperature or 
emission measure with kernel size or altitude. In this case plasma density 
calculated from A^e = \/EM/V should systematically decrease with kernel 
altitude because of the increasing kernel volume. The same relation should 
be observed for plasma pressure because the temperature in our sample is 
also uncorrelated with geometrical aspects of flare kernels. As mentioned 
above, the kernel volume rises with altitude as F oc h}'^^ , which should 
result in density or pressure drop with altitude as 

AT, = {EM/Vfl'' oc /i-0-85±o.03 

and 

Pe oc iVe oc h'^'^^ . 

This slope is consistent with the observed relation in the sub-sample of M- 
class flares. For C and X-class flares the agreement is worse, but this is 
explained by the X-class sub-sample having small number of events, and the 
sub-sample of C-class kernel being ill-spaced as mentioned above. 

We interpret these relations as nearly identical concerning the slope, but 
differing in the intercept parameter. Stronger flares involve more emission 
measure within the same volume, i.e. larger flares need higher densities and 
pressures. 



5. Energy release within the kernels 

Loop-top kernels are the brightest parts of the soft X-ray flares and pre- 
sumably are the places where the plasma heating is most effective. The 
next physical parameter we estimated for the flare kernels was the heating 
rate during the kernel maximum. Close to flare maximum, when the energy 
transport by plasma flow is negligible we can build a simplified model of the 
evolution of thermal energy {Eth) content as the sum of heating {Eh) and 
cooling processes (radiative and conductive losses: Er and Eq)'- 

-jrEth = Eh — Ec — Er . 
at 

Moreover, at the flare maximum the changes of thermal energy content 
are insignificant, so we can estimate the energy release rate for this moment 
as 
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Eff ~ Ec + Er . 

Radiative loss can be calculated knowing the temperature and plasma 
density estimates, Er = N^A{T), where A(T) is the radiative loss function 
taken e.g. from Rosner, Tucker and Vaiana (1978). 

Wc estimate the kernel heating in the way already shown by Jakimiec, 
Falewicz and Tomczak (2002), which we briefly summarize here. We assume 
that the kernel is nearly spherical and it is heated by some process with the 
mean rate Eff. The total flux of thermal energy flowing out of the kernel may 
be estimated as (47r/3) i?^(i?iy — Er), where R is the kernel radius and Er is 
the mean radiative loss from the kernel. This flux has to be transported down 
by thermal conduction along the both legs of the loop containing the kernel. 
If we take the diameter of this loop to be the same as that of the kernel, 
then the conductive flux in the legs is Fq = {2/3)R{Eh — Er). Knowing 
this value we can integrate the equation of conduction, 

Fc = KoT'^-^dt/ds, 
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Figure 8. Relation between the plasma pressure in the kernel and its altitude above the 
photosphere. The lines show power-law fits to each of three GOES sub-samples. 



along the loop of semilength L, assuming constant Fc- The result is 

0.286r3-5 = {Fc/ko)L, 

where T is now the temperature of the kernel. Prom these two equations we 
can estimate the mean heating within the kernel as 

Eh = 3.9 X 10-'^ T^-^/RL + Er. 

The first term on the right hand of this formula is the conductive loss from 
the kernel. Its formulation differs than the radiative losses from the simple 
loop. We want to stress that this estimation does not implicitly involve which 
process is responsible for the kernel formation. 

We applied the above formula to the analyzed sample, taking the geo- 
metrical mean of two semi-axes seen on SXT image as the kernel radius, and 
the approximation of the loop semi-length as L = {n/2) h. A typical value 
of heating is a few ergs per second in each cubic centimeter; however, the 
observed values within our sample span two orders of magnitude. Figure 9 
shows a clear decrease of the heating rate with the altitude of a kernel. 
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C and M-class kernels show a very similar dependence of heating with 
kernel altitude, log Eh = 15.63 — (1.70 it 0.05) log/i. This is similar to the 
results obtained by Jakimiec, Falewicz and Tomczak (2002) who analyzed 
a set of 27 limb flares of mostly M-class and loop-length longer than 10^ 
cm. They estimated that the flare heating rate drops with loop-length as 
log Eh = 13.91 — 1.46 log L. Bcik-St§slicka and Jakimiec (2005) enlarged this 
sample by adding 10 slow long duration flares, which allowed them to widen 
the range of analyzed loop-lengths by an order of magnitude. They obtained 
the relation log Eh = 16.3 — 1.72 log L. Within errors this is nearly identical 
to the relation obtained here for C and M-class flares. Our sample consists 
also of 8 X-class flares, for which the dependence between heating rate and 
kernel altitude is log Eh = 13.0 — (1.33 it 0.20) logh, but the same slope 
as in C and M flares cannot be excluded. X-class flares require about four 
times higher heating rates, which in our calculations is caused mainly by 
distinct increase of emission measure. The radiative losses in these kernels 
are correspondingly high to balance the higher level of heating. 

These slopes again seem to result mainly from the geometrical properties 
of the kernels. During the flare maximum radiative and conductive losses 
are usually comparable, Ec ~ Er (see Fig. 10). With the lack of any 
dependence of temperature on kernel altitude, the conductive losses scale as 
Ec ~ {RL)-^. Taking our result, R ~ L°-^^ we obtain Eh Ec L'^-^^. 
This slope differs from the relations in Fig. 9 by no more than 3(j. 

Figure 10 also shows that the ratio between the radiative and conduc- 
tive losses changes when going from weak to strong flares. The stronger 
the flare, the more important the influence of radiative losses. Conduc- 
tive losses depend on the geometrical aspects and the plasma temperature. 
The geometrical properties of flare kernels arc more or less the same for 
each GOES sub-sample. Taking the relation between temperature and emis- 
sion measure observed in our sample we can show that Ec oc T'^-^/RL oc 
EM^-^^^ /h^-^^. Radiative losses are proportional to emission measure, Er = 
EM X A{T)/V oc EM^-^^/h^-^^, where we assumed that A(r) oc T'^/s. 
Both losses have similar dependence on the kernel altitTide, but diflFerent 
on the emission measure, which makes radiative cooling more important in 
strong flares, Er/Ec oc EM^-^^ /hP'^^ . This result is, however, sensitive to 
the assumed dependence between temperature and emission measure. The 
SXT instrument has limited sensitivity to the hot plasma which is the cause 
of the quite flat relation between these two parameters in our sample. It 
is enough to assume relation T oc EM^''^'^, which is power-law equivalent 
of the relation presented by Feldman et al. (1996), to get the ratio Er/Ec 
independent of the kernel emission measure. 

The limited temperature sensitivity of the SXT affects both radiative 
and conductive cooling. Underestimation of kernel temperature results in 
overestimation of radiative and underestimation of conductive losses. The 
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second effect is stronger, because it depends on temperature to a much 
higher power. Thus the heating rate calculated in this paper seems also to 
be underestimated and should be treated as a lower limit. This bias should be 
more important for stronger flares. A typical value of kernel temperatures for 
X-class flares in our sample is T ~ 12 MK, which seems to be a factor 2 less 
than temperatures reported by other similar instruments. This would result 
in underestimating of conductive losses by about one order of magnitude 
and overestimating the radiative losses by a factor 1.6. Good estimation 
of cooling and heating rates for strong flares requires an imaging X-ray 
telescope with better temperature sensitivity than SXT. 



6. Discussion and conclusions 

The observed relation between kernel size and its altitude is not an obvious 
result. If we assume that kernel is part of the loop that appears to contain 
it, than the loop cross-section area must show the same relation. This means 
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Figure 10. Comparison of the radiative and conductive losses within the analyzed flare 
kernels. The solid line shows equality between the two quantities, while the two dashed 
lines mark the borders when one loss exceeds the other by a factor 10. 



that the higher flare loops arc not hnearly scaled version of the smaller loops, 
i.e. the loop diameter, $ w (4A/7r)^/^ is not proportional to the loop length, 
but rather to its square root, $ ~ /ji.i3/2 ^ ^0.56 fpj^g higher loops must be 
systematically more narrow than the smaller ones. 

The relation between plasma pressure and the kernel altitude allows us 
to put some constrains on the magnetic field structure in active regions. 
The kernel plasma of a given pressure must be trapped within a field with 
magnetic pressure greater than plasma pressure, B'^/Stt > pe- This allows 
us to determine the minimum magnetic strength necessary to contain the 
kernel plasma as log = 5.86 — 0.43 log h for M-class kernels. This con- 
straint slightly depends on the GOES class. If wc assume that for all flares 
logpe decreases with altitude with the same slope as for M-class kernels, 
the intercept parameter of log Bj^in in C-class flares is 5.71 while in X-class 
flares it is 6.11. This estimation does not take into account the influence of 
possible field helicity, which allows a weaker field effectively to trap a kernel 
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Figure 11. Relation between the magnetic field necessary to trap the kernel plasma and 
the kernel altitude. The two Unes presents the dipole field, according to Aschwanden et al. 
(1999) approximation with ho = 7.5 x 10® cm and Bo equal 70 and 400 G. 



plasma. The field helicity is unfortunately difficult to determine with present 
X-ray observations of flares. 

The observed relation Bmin{h) is also in quite good agreement with 
estimates of magnetic field in active regions shown by Aschwanden et al. 
(1999). They approximated the field with a magnetic dipole, B{h) ~ So(l-|- 
h/hf))^'^, where hf) is the dipole depth estimated by them as 7.5 x 10^ cm, 
and Bq is the mean magnetic field at the photospheric level. Fig. 11 shows 
that our sample is consistent with the dipole field with Bq in the range of 
70 to 400 G. 

The relation obtained between kernel heating and its altitude may also 
shed light on the magnetic field structure in active regions. If we hypotheti- 
cally assume, as in Jakimiec, Falewicz and Tomczak (2002), that the energy 
release rate is somehow proportional to the density of magnetic energy 
Eh ~ B^/Sn, then we should observe decrease of the magnetic field with 
height as 5 ~ Jakimiec, Falewicz and Tomczak (2002) suggested in 

the same way that magnetic field responsible for fiare heating decreases with 
height as B ^ This relation has a different slope from the decrease 
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of magnetic field necessary to trap the kernel plasma, but both relations 
may describe different parts of the field supporting the kernel. The field 
necessary to contain the kernel plasma may be more related to the kernel 
surroundings, while the field responsible for heating is obviously related to 
the location of energy release. These two places are not by definition the 
same location as it is discussed e.g. in Karlicky and Barta (2006). 

The reported here relation between the kernel size and its altitude should 
be considered regarding the flare scaling laws. It is generally accepted that 
during the flare maximum radiative and conductive losses are of the same 
order, ~ Ec- Conductive losses for a simple flaring loop without a kernel, 
Ec oc T^'^/L^, changes into Eq oc T^'^/RL when we describe a loop-top 
kernel. In this case the scaling takes the form iV|A(T) ~ T^-^/RL. If we 
substitute A(r) ~ T~/^, then we can rewrite this scaling as 

The relation between kernel size and its altitude is i? ~ /i^'^^, which makes 
the product RL ^ h^'^^ close to kernel volume R^ ^ h^-^^. Thus, the left side 
of the above equation roughly scales as kernel emission measure. The scaling 
Er ~ Ec turns into a relation between flare temperature and emission 
measure, EM ~ T^-^+P. Taking /? = 2/3 we achieve T ~ EM^-"^"^ which is 
very close to the relation shown by Feldman et al. (1996). 

The existence of flare kernels should be taken into account also when 
analyzing stellar flares. It is hard to imagine that in stellar flares loop-top 
kernels arc not present when in solar ones they arc common. We should be 
careful when transferring the scaling laws of solar kernels onto flares on other 
stars. Recognizing how the overall magnetic structure affects formation and 
evolution of a flare kernel is necessary. The Sun represents stars with very 
low magnetic activity. We cannot exclude the possibility that the kernel 
scaling laws may have different formulation on highly active stars, which 
seems to be fully covered by active regions. However, we remark that the 
similarity of a T — EM relation for stellar flares with a solar one (see e.g. 
Fig. 1 in Shibata and Yokoyama, 2002) suggests that kernel scaling laws on 
other stars may not differ substantially. 

Flare loop-top kernels are obvious characteristics of SXR solar flares. 
The kernel emission is usually constant part of the total flare emission, 
what allows to describe overall evolution of a flare by the evolution of its 
kernel(s). A correct model of either solar or stellar flares cannot be built 
without understanding their properties and nature. 

This study shows that the set of basic obscrvablcs describing the flare 
loop-top kernels, which we obtain from SXT images, separates into two inde- 
pendent pairs. One pair describes the geometrical properties of a kernel (size 
and altitude), the other pair is plasma temperature and emission measure. 
As we report here, kernel size and altitude are not independent parameters. 
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Although in individual flares one can see different types of relation between 
them, when comparing many flares, a clear dependence emerges: the higher 
the kernel, the bigger the size. The second pair, temperature and emission 
measure, is also internally related, but these parameters are not correlated 
with the kernel altitude and size. To describe many flares, it is sufficient to 
take only two observables, viz. emission measure and kernel height. Other 
physical parameters of emitting plasma like density, pressure or heating rate 
are related to these two quantities. 

To determine more precisely the kernel scaling laws, we need a larger 
sample of limb flares. We expect that the forthcoming imaging instrument 
in soft X-rays, XRT onboard Solar-B, will open new possibilities in ana- 
lyzing the physics of flare loop-top kernels. The expected improved angular 
resolution will allow better determination of kernel geometrical properties. 
Improved thermal sensitivity of XRT shoTild allow us to achieve more precise 
determination of physical properties of the emitting plasma. We may expect 
better determination of conductive losses and heating rates, which strongly 
depend on the estimated kernel temperature. The next few years of soft X- 
ray observation of the Sun will allow us to substantially enlarge the sample 
of suitable, near-limb flares and take a closer look on the kernel properties. 
Difficult to model, flare kernels still await our attention. 
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